Reflection of short polarized optical pulses from periodic and aperiodic 

multiple quantum well structures 
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We study the reflection of polarized optical pulses from resonant photonic structures formed by 
periodic, Fibonacci, and gradient sequences of quantum wells. The form and polarization of the 
reflected pulse are shown to be determined by the structure design and optical length. In structures 
with periodic quantum well arrangement, the response to ultrashort pulse is an optical signal with a 
sharp rise followed by an exponential decay or Bessel beats depending on the structure length. The 
duration of reflected pulses non-monotonically depends on the number of quantum wells reaching 
the minimum for a certain structure length which corresponds to the transition from superradiant to 
photonic-crystalline regime. We also study the conversion of pulse polarization in the longitudinal 
external magnetic field which splits the exciton resonance. Comparing periodic, Fibonacci, and 
gradient structures we show that the latter are more efficient for the conversion from linear to 
circular polarization. 

PACS numbers: 42.70.Qs, 78.47.jg, 78.67.Dc, 71.35.-y 
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I. INTRODUCTION 

The optical spectroscopy of photonic structures, where 
the susceptibility is spatially modulated on the scale 
of light wavelength, has been attracting a lot of atten- 
tion due to fundamental interest in the light-matter cou- 
pling as well as applications of such systems to opti- 
cal devices J^- The photonic structures sharing similar 
physics can be realized in different systems ranging from 
crystals for 7-rays,^ solid state heterostructures for in- 
frared, visible, or ultraviolet ligln— £ to the optical lattices 
of laser-cooled atoms 7 - and the chains of optical cavities^ 
Of special interest in this field are the resonant photonic 
structures based on periodically arranged semiconductor 
quantum wells (QWs) or dots (QDs) where the light cou- 
ples to excitons confined in the QWs or QDs leading to 
formation of polaritonic modes.— An advantage of such 
structures is the possibility of electric control of the light- 
exciton coupling as well as tuning the spectral position 
of the exciton resonance . 10 ' 11 

The QW-based resonant photonic crystals suggested 
two decades ago^ are extensively studied at present both 
theoretically^— and experimentall y 1 1 1 18 " — with the fo- 
cus on their spectral characteristics. The progress in the 
optical spectroscopy technique makes it possible to study 
the response of photonic structures at the pico- and fem- 
tosecond time scale ; 24 ' 25 which stimulates theoretical in- 
vestigations of ultrafast optical processes. In the present 
paper, we develop the theory of time- and polarization- 
resolved linear response of multiple QW systems. We 
show that the optical spectroscopy in the time domain, 
being complementary to that in the frequency domain, 
provides direct information on the lifetime and coher- 
ence of optical excitations, which cannot be easily ob- 
tained from the stationary reflection spectra. Previous 
studies of time-dependent reflectivity were limited to pe- 
riodic QW structures i 13 ' 24 ' 26 Here, we focus on the effect 
of geometrical arrangement of QWs on the optical re- 



sponse and compare the responses of periodic, quasipcri- 
odic Fibonacc i 27 ' 28 as well as gradient QW structures. 
We also analyze the linear-to-linear and linear-to-circular 
polarization conversion of pulses which occurs in the ex- 
ternal magnetic field. 

The paper is organized as follows. The multiple QW 
structures and the general expressions for the structure 
optical response are presented in Sec. [TTJ In Sees. IIIII 
and IIV| we study the reflection of short optical pulses 
from periodic and aperiodic Bragg structures, respec- 
tively. The polarization conversion is analyzed in Sec.fVl 
The results of the paper are summarized in Sec. IVII 



II. MODEL 

We consider the layered heterostructure which consists 
of N identical quantum wells embedded in a dielectric 
matrix, see Fig. [TJ The wells are centered at the points 
Zj (j = 1 . . . N) with z\ = 0. The structure is charac- 
terized by the background refractive index rib- Addition- 
ally, the linear response of each QW to the electromag- 
netic field has the exciton resonance at the frequency ujq. 
The corresponding frequency dependence of the ampli- 
tude reflection and transmission coefficients, r\(oS) and 
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FIG. 1. (Color online) Illustration of optical pulse reflection 
from a multiple quantum well structure. 
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ti(u>), respectively, is given by 



n(w) = - 



ir r 



u - uj + i(r + r) 



ti(w) = l + ri(w), (1) 



where To and T are the radiative and non-radiative de- 
cay rates. We assume that the QWs are all identical 
and uniform in the interface plane, so that the effects of 
inhomogeneity2£ can be neglected. 

The light pulse is normally incident on the QW struc- 
ture along the z axis and reflected. In the left half-space, 
z < 0, the electric field E(z, t) is the sum of incident and 
reflected waves, 

E(z, t) = Eo(t- zn b /c) + E R (t + zn b /c) , (2) 

where c is the speed of light. The reflected pulse Er{£) 
is determined by the linear response function p{r) as fol- 
lows: 



/•OO 

E R {t) = / P (T)E (t-T)dT . 

Jo 



(3) 



The function p(r), in turn, represents the Fourier trans- 
formation of the reflection coefficient of the structure in 
the frequency domain r(w), 



P(t) 



r(u)e~ iut du/(2ir) . 



(4) 



We calculate the reflected pulse En{t) after Eqs. §3§ 
and ((!]). The reflection coefficient for an arbitrary QW 
structure is found by the transfer-matrix technique^ In 
this approach, each layer of the structure is described by 
the transfer matrix which connects the amplitudes of in- 
coming and outgoing waves. The corresponding matrices 
for an individual QW and interwell distance d are given 
by 



Tqw 



ti(w) 



t?H-r?H n(w) 
-n(w) 1 



cxp(iujdn b /c) 

exp(— iLodrifj/c) 



(5) 



(6) 



respectively. The total transfer matrix of the structure T 
is given by the product of the individual transfer matrices 
of all layers. It enables one to calculate the reflection 
r{u>) and transmission t(uj) coefficients from the matrix 
equation 



T 



1 




t(w)" 


r(uj)_ 








(7) 



In what follows we are interested in the optical re- 
sponse to ultrashort pulses which have the carrier fre- 
quency at the resonance frequency wo and the spectral 
width 5uj much broader than the scale of the exciton- 
rclatcd spectral features, <C Suj <C loq. Since the 

typical band width of exciton reflection in periodic QW 



structures is of the order of few meV ) 11 ^ 4 ' 28 this con- 
dition can be realized for sub-picosecond pulses. Then, 
the time dependence of the reflected pulse electric field 
has the form E R (t) = p(t) J+°° E { T )e luloT dr and follows 
p(t). 



III. PERIODIC BRAGG STRUCTURES 

In this Section, we consider periodic QW structures 
where the interwell distance d equals to half of the light 
wavelength at the exciton resonance frequency 



irc/(ui n b ) 



(8) 



The calculation by the transfer matrix technique shows 
that the reflection coefficient of such structure containing 
N QWs in the frequency range \u> — <C loq can be 
compactly presented by^^ 



rjv(w) = ( 2ww 



lcot<I> 



(9) 



Here, w = (u-uj )/A, w = w + iT/A, A = ^(2/ir)uj r a 
is the halfwidth of the stop band, $ = Iv^J w 2 — w/w 
is a resonant contribution to light phase incursion, and 
v = NTq/ A is the normalized number of QWs. 

Figure [2] shows the envelopes of reflected pulses in 
the time domain calculated numerically for structures 
with the various QW number. The curves in Figs. [Ha) 
and [2]Jb) are plotted, respectively, for the cases of zero 
and rather strong non-radiative decay. One can see that 
the shapes of reflected signals in both cases are similar 
and drastically depend on the QW structure length. In 
short structures, the pulse envelope has a sharp rise fol- 
lowed by a slow monotonous decay. Both the signal am- 
plitude and the decay rate increase with the QW num- 
ber N. In contrast, the reflected pulse envelope in long 
and, particular, semi-infinite structures exhibits oscilla- 
tions. The period of the oscillations is determined by the 
stop band width A. The calculation shows that, in the 
absence of non-radiative decay, the transition between 
these regimes occurs at v sa 0.66, which corresponds to 
N = 106 for QW parameters given in the caption of 
Fig-HHa). For comparison, we plot the reflection spectra 
|r(o;)| 2 of the QW structures in the insets of Figs. EJa) 
and[2]Jb). With the increase in QW number, the reflec- 
tion spectrum gradually evolves from the Lorentzian to 
a silk-hat profile without prominent changes. It demon- 
strates that the time-resolved spectroscopy can be an ef- 
ficient tool to study photonic structures supplementing 
the frequency-resolved approach. In the rest of Section, 
we present analytical results for short and long periodic 
structures and study the transition between these two 
regimes of light reflection. We also analyze the role of 
non-radiative decay. 



3 



infinite 




FIG. 2. (Color online) The envelopes of reflected pulses calcu- 
lated for the resonant Bragg structures with various numbers 
of QWs and cjo/Fo = 4 x 10 4 . Figures (a) and (b) corre- 
spond to the absence, F = 0, and rather strong, V — 20Fo, 
exciton non-radiative decay, respectively. The dashed curves 
depict the reflected pulse envelopes for the semi-infinite Bragg 
structure. Insets show the reflection spectra of the same QW 
structures in the frequency domain. 



Transition from superradiant to 
photonic-crystalline regime 



The optical response of resonant Bragg structures with 
the small number of layers, v ^ 1 or TV <C ^ljq/Tq ~ 100 
for typical QWs, is described by two times: the short 
flight time of light through the structure and back t[ = 
2irN/uio and the long lifetime of QW excitation. Dur- 
ing its flight through the structure, the incident pulse 
sequentially excites QWs. Therefore, the reflected sig- 
nal amplitude linearly increases with time at i < if, see 
Figs. [2ja) and [2{b). The signal behavior at the times 
t 3> tf is determined by the excitation decay. To describe 
it analytically one should consider the limit \w\, \ w\ <C 1 
in Eq. which yields the Lorcntzian for the reflection 



coefficient in the frequency domain^ 

iNTp 

r N {uj) = , ., Mr , -ps • (10) 

lj - uq + i(at + r) 

The corresponding response function /9jy(i) determining 
the time dependence of reflected signal has the form 

PN (t) = -AToe-^+^'e- 1 " * 0(t) , (11) 

where 9(t) is the Heaviside step function. The pulse en- 
velope exponentially decays with time at the rate pro- 
portional to the QW numbcri 13 ' 26 Such an enhancement 
of the decay rate is a direct manifestation of the Dicke 
supcrradiancei 29 ' 30 

The further growth of the structure length leads to 
the increase of the flight time and decrease of the decay 
time. At N ~ ycJo/Trj, both times become compara- 
ble to each other and the optical response qualitatively 
changes. This corresponds to the transition from super- 
radiant to photonic-crystalline regime with the stop band 
formation. In the limiting case of N ^> ^Jujq/Tq, the 
structure reflectivity |r(uj)| 2 is close to unity in the stop 
band spectral range \w — u>o\ < A and rapidly decreases 
outside this stop band, see inset in Fig. [U[a). The re- 
flection coefficient of the semi-infinite structure is given 
by 



?"oo (w) 



1 



(\/wW + yjww — 1^ 



(12) 



The corresponding response function in the time domain 
at T < A and t > T 2 /A 3 has the form 



Poo(t) 



J 2 (At) 



^(rt/2) 



-rt/2 e -iu o t0( t ) ; 



(13) 

where J 2 and I\ are the Bessel and modified Bessel func- 
tions of the first kind. In particular, in the absence of 
non-radiative decay, (t) is simplified to 



p 00 (t) = --MtA)e- i "° t 6(t) 



t 



(14) 



Equations (| 13[) and (fT4"l) demonstrate that the envelope of 
the pulse reflected from a long structure oscillates in time. 
The oscillation frequency is determined by the stop band 
halfwidth A. We note that oscillations of such kind are 
known for the scattering of 7-rays by nuclei^ However, to 
the best of our knowledge, they have not been considered 
so far for multiple QW structures. 

The transition from superradiant to photonic- 
crystalline regime can be also revealed by studying the 
normalized energy Vn and duration Tn of the reflected 
pulse which we define by 

/>OC 1 />oo 

V N = \ P (t)\ 2 dt, T N = — \ P (t)\ 2 tdt. (15) 

Jo l J N Jo 

The dependences of Vn and T/v on the number of QWs N 
for periodic structures are shown in Figs.^a) and[3Jb), 
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100 200 

Number of QWs N 

FIG. 3. (Color online) Dependences of the normalized energy 
Vn [figure (a)] and duration Tn [figure (b)] of the reflected 
pulse on the number of QWs. The curves are plotted for the 
periodic Bragg structures and the Fibonacci structures with 
different Bragg conditions, wo/ T = 4 x 10 4 , and T = 20F . 



respectively, by black curves. The reflected energy in- 
creases with the QW number in the superradiant regime, 
V N = ^V 2 rg/[2(iVr + T)], and reaches the value of 
Voo = 16A/(157r) in the photonic-crystalline regime pro- 
vided r <C A. The transition from the linear growth 
to saturation of Vn corresponds to the transition from 
superradiant to photonic-crystalline regime. The depen- 
dence of the reflected pulse duration Tn on the QW 
number N is even more pronounced. In the superra- 
diant regime, it decreases with QW number as Tn = 
l/(2ATo + 2r). In the area of the superradiant-photonic- 
crystallinc regime transition Tn reaches the minimum 
sa 2.1/A at N sa 0.51A/r (for r < A) and then 
again increases with the structure length. Finally, it sat- 
urates at the value Too = 157r/(16A) (for T <C A) for 
the semi-infinite structure. Numerical calculation shows 
that, for the QW parameters given in the caption to 
Fig. [3l the shortest time Tn is achieved in the struc- 
ture with N = 58. This provides a receipt for designing 
the fast-reflecting Bragg structures: the structure length 
should be intermediate and correspond to the transition 
between superradiant and photonic-crystalline regimes. 



B. Long structures. Echo and Borrmann effects. 

Now we discuss the features of pulse reflection from 
long, N ^> ^/wo/r , but finite QW structures. The re- 
flection coefficient of such structures rjv(w) is described 
by Eq. (|12l) with corrections which can be obtained by 
expanding Eq. Q in the geometrical series in the param- 



eter exp(2ii>). Such a procedure yields 

oo 

p=i 



(16) 



where 



4V WWy/wW — 1 



'm + s/WW — 1 



4p 



cxp(2i$p) . (17) 



The corresponding response function in the time domain 
Piv(i) has the form 



P =i 



(18) 



with 5p${t) being the Fourier image of 6r%>(uj). The 

correction 5pN (t) can be interpreted as the contribu- 
tion to the response function caused by p-th reflection 
of the optical pulse from the structure back edge. In- 
deed, 8p^\t) 7^ only for t > ptf, where ptf = 2ttNp/uj 
is the time required for the light to travel through the 
structure and back p times. In particular, in the absence 
of non-radiative decay, the corrections 5p N (t) have the 
approximate form 



(p), 



Xt-pttfy-y^pttft 
' (t + puyp+v* 



Ay^ 2 - (pi f ) : 




Ptf), (19) 



which is obtained by the stationary phase method. 

The first-order reflection from the structure back edge 
is clearly seen in the pulse envelopes plotted for the QW 
structures with N = 300 in Figs. Ufa) and^b). It ap- 
pears at the time t = tf [tA w 7.5 in Figs.[2fa) and[2fb)] 
leading to an uprise of the pulse envelope. This echo-like 
feature is more pronounced in Fig. Ufa) where we com- 
pare the response functions /Ojv(i) for long (N = 10 3 or 
v Ri 6.3, red solid curve) and semi-infinite (black dashed 
curve) QW structures. The curves in Fig. Ufa) are all 
calculated for structures with the strong non-radiative 
decay. 

Another interesting difference between the response 
functions of a long finite and the semi-infinite QW struc- 
tures is their behavior at large times, compare solid and 
dashed curves in Fig. Ufa) . The response function of the 
semi-infinite structure in the presence of non-radiative 
decay is described by Eq. (|13|) . The counterintuitive fea- 
ture of Eq. ([13]) is a power decay of poa (t) at large times 
that is clearly seen in the double logarithmic scale in 
the inset of Fig. Ufa) . Instead of the exponential decay, 
expected for lossy systems, the time response function 
has the asymptotics Poo(t) w \/T / (7rA 2 t 3 ) at t ^> 1/r. 
Such a behavior is closely related to the Borrmann effect, 
known in the X-ray physics as the origin of extraordinary 
transmission of radiation through absorptive mcdiai 3 ^ In 
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FIG. 4. (Color online) (a) The envelopes of the pulses re- 
flected from the periodic Bragg structure with 1000 QWs (red 
solid curve) and semi- infinite structure (black dashed curve). 
The curves are calculated for io /T = 4 x 10 4 and T = 20r . 
The inset shows the envelopes plotted in the double logarith- 
mic scale, (b) Reflection spectra of the same QW structures 
in the frequency domain. The inset shows the reflection spec- 
trum details at ui ~ cjq. 



our case of resonant Bragg structures, the effect comes 
from the suppression of radiation absorption at the res- 
onant frequency ujq. Indeed, in the semi- infinite struc- 
ture, the steady-state distribution of radiation electric 
field at the resonance frequency is a standing wave with 
the nodes situated at the QW centers^ Such a wave 
does not interact with excitons in QWs. This leads to 
a narrow dip in the absorption spectrum and a narrow 
peak in the reflection spectrum |r(w)| 2 at the frequency 
wo with the characteristic width of T^ see Fig.@|b) and 
the inset. In the semi-infinite QW structure, |r(wo)| = 1, 
which results in a power decay of Poo(t) at large times. 

In a long but finite structure, the peak in the reflec- 
tion spectrum remains, however, |rjv(w)| does not reach 
unity, see the inset in Fig.|4fb). Therefore, a part of the 
response function with a power-like decay can be still dis- 
tinguished but the behavior of pNit) at t — >• oo is always 
described by the exponent. Such a response function is 



shown in Fig. QJa) by red solid curve. The decay rate of 
PN(t) is determined by the pole of Eq. © with the small- 
est negative imaginary part, which yields the asymptotics 
p N {t) w 7r 2 r 2 /(2i' 3 A) exp(-7r 2 rt/2i/) at t > v^/T. 



IV. APERIODIC BRAGG STRUCTURES 

In this section, we analyze the effects of non-periodicity 
of QW arrangement and detuning from the Bragg condi- 
tion on the optical response. We focus on quasiperiodic 
Fibonacci structures (Sec. IIV Ap and gradient structures 
where the interwell distance gradually changes from the 
structure forefront to its back edge (Sec. HVBl) . 



A. Fibonacci structures 

In canonic Fibonacci QW chains, the interwell dis- 
tances can take one of the two values, a or b, with a/b 
being equal to the golden ratio r = (a/5 + l)/2 « i.62. 
The sequence of the distances is determined by the re- 
current relation similar to that for the Fibonacci num- 
bers. The sequence of the order m, J- m , is given by 
J r m _iJ r m _2 with T\ = b and Ti = a.— In particular, 
the chain = abaab describes the structure consisting 
of N = 6 quantum wells positioned at z\ — 0, Z2 = a, 
23 = a + b, Z4 = 2a + 6, 2:5 = 3a + 6, and zq = 3a + 26. 

Despite the non-periodicity, the Fibonacci QW struc- 
tures can be tuned to the resonant reflection.— The cor- 
responding Bragg condition for such structures has the 
form 



w n b - I ti 
-a = 7r h H 



(20) 



where d = (3 — r)b is the mean interwell distance and h 
and h! are the integer numbers. In contrast to periodic 
systems, the Bragg resonances in the Fibonacci struc- 
tures are characterized by two numbers h and h' , which 
reflects the fact that these structures can be obtained 
by the cut-and-project method from the two-dimensional 
square lattice! 33 ' 34 

The envelopes of reflected pulses in the time domain 
for the Fibonacci structures tuned to the Bragg reso- 
nance (h, h') = (1, 0) are shown in Fig. [5] by solid curves. 
Figures [5ja) and EJb) correspond to the cases of weak 
and strong non-radiative decay, respectively. The curves 
are calculated numerically by the transfer matrix tech- 
nique. One can see that the reflected pulses exhibit os- 
cillations even for structures with rather small QW num- 
bers N which operate in the superradiant regime. This 
is in contrast to periodic Bragg structures where the os- 
cillations emerge in the photonic-crystalline regime, i.e., 
at large N, see Fig. [2j The origin of such a behavior is 
a dip in the reflection spectrum of the Fibonacci struc- 
tures at lu w ljo, see insets in Figs. [SJa) and[SJb), which 
is caused by non-periodicity of the QW arrangement i22 
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FIG. 5. (Color online) The envelopes of reflected pulses cal- 
culated for the Fibonacci structures with various numbers 
of QWs, the Bragg condition (h,h') = (1,0), and luo/Tq — 
4x 10 4 . Fi gures (a) and (b) correspond to weak, F = Fo, and 
strong, r = 20ro, non-radiative decay, respectively. Black 
dashed curves depict the pulse envelopes for the semi-infinite 
Fibonacci structure. Black dotted curves present the analyti- 
cal asymptotics Eq. ()25[) . Red dash-dotted curves are plotted 
after analytical Eq. (|24l) obtained for the superradiant regime. 
The insets show the reflection spectra of the same QW struc- 
tures in the frequency domain. 



Similar spectrum structure has been observed for peri- 
odic Bragg systems with two QWs in the unit cell^ and 
for layered nuclear systems^ 

To analyze the reflection of optical pulses from the Fi- 
bonacci structures in more detail we use the two-wave 
approximation. In this approach, the eigen electromag- 
netic states of the structure in the vicinity of the Bragg 
resonance (h, h') are approximated by the sum of two 
plane waves with the wave vectors K and K — Ghh' , with 
Ghh' = (2ir/d){h + h'/r) being the diffraction vector^ 
The approach yields the reflection spectra of the Fi- 
bonacci structures 



1 1 h' 



UJq - LO - X 



if2cot(2O-/Vr /A£ fe ,) ' 



(21) 



where 



X = 



A 2 

^hh' 



2(u}q — ui — iT) 



n = V(x + u-u ) 2 -x 2 \h h '\ 2 , 



Ahh' = A/ y/h + hi Jt, and fhh' is the structure factor- 
determining the strength of the Bragg diffraction peak, 



fhh' 



sin Shh > 
Shh' 



■ exp 



■S 



TTrirh' — h) 



T 2 + 1 ' 

(22) 

In a simple periodic structure fhh' — 1 an d Eq. (|2Ijl re- 
duces to Eq. ([9]), while in quasicrystals \fhh'\ < F Com- 
parison of the reflection spectra given by Eq. (f2~Tj) with 
those obtained numerically by the transfer matrix tech- 
nique shows that Eq. (|21[) is valid in the whole frequency 
range of the resonant reflection except for a narrow cen- 
tral region. 

In the regime of superradiancc, Eq. (|21[) assumes the 
form 



r N (uj) 



hh' 



1 + \OLhh' cot 



ahh' NT o 
luq — ui — iT 



(23) 



where ahh' = \/5 — |A/i'| 2 characterizes the strength 
of aperiodicity. For a periodic structure ahh' = and 
Eq. p3|) reduces to Eq. (fl"0|). The reflection spectrum 
given by Eq. (|23|) has a wide peak at the frequency wo 
with the width proportional to NTq, which is a signa- 
ture of the superradiant regime. In the middle of the 
peak, there is a structured dip of the width ~ ahh' NTq, 
which is caused by the spacial nonperiodicity of the sys- 
tem. Such a multipeak structure of the reflection spec- 
trum rjv(w) leads, in turn, to a beating pattern in the 
response function pjv(i) in the time domain, see also 
Ref. [37l . To obtain analytical expression for pj^(t), we 
expand Eq. (|23[) in the Taylor series at ahh' ~^ 1 an d 
perform the Fourier transformation Eq. (U). This proce- 
dure yields 



p N (t) = -2f hh 'J2 



fe-i 



k=l 



OkhA^ - ahh') 
{l + a hh ') k+1 



2ka hh 'NT 



t 



x Jr( V8^WVIV) e- rt e- iuot 9(t) . (24) 



The response functions given by Eq. (|24[) for the QW 
number N ~ 60 and different non-radiative decay rates 
are plotted by dash-dotted curves in Figs.[5ja) and[5jb). 
One can see that the equation describes the exact solu- 
tions shown by solid curves quite well. We note that for 
small numbers of h and hi , where the parameter ahh' is 
close to unity, the dominant contribution to pn (t) comes 
from the first term (k = 1) of the series (|24[) . 

In the opposite case of very long structures, N 
y/uj o /T , the multiple reflection of light from quantum 
wells in the Fibonacci structure leads to the formation of 
two wide stop bands, see insets in Fig.[5ja) andJSJb). The 



corresponding reflection coefficient roo(w) for the semi- 
infinite structure is obtained from Eq. (f2"Tj) by the re- 
placement of cot^ilATo/A^,) in the denominator with 
(— i)sign(InrO). The subsequent Fourier transformation 
of Too (ui) by the stationary phase method yields the re- 
sponse function in the time domain 



P°o(t) « // 



hh' 



\ 3 K, h >t 3 



„-rt/2 - 
c ' e 



6{t) 















X 


^/cj + sin (uj+t - 


-j) 




\J(jJ- cos (ui-t - 


-;)] 



which is valid at t > 



(25) 



Here, the frequencies 



w ± = JA 2 h h ,{l ± |/wl)/ 2 - T 2 / 4 determine the inner 
and outer edges of the stop bands measured from tu . The 
analytical function (|25|) shown in Figs. [5{a) and 03b) by 
dotted curves is in a perfect agreement with the results 
of numerical calculation by the transfer matrix technique 
(dashed curves) at large times. 

The transition from superradiant to photonic- 
crystalline regime with the increase of the QW number 
leads to a change in the normalized energy Vn and du- 
ration Tn 01 the reflected pulse, Eq. ([To]). The depen- 
dences of Vn and Tn on the number of QWs for the Fi- 
bonacci structures with different indices (h, h') are shown 
in Fig.^a) and[3jb), respectively They are qualitatively 
similar to those for periodic structures. The energy Vn 
increases linearly with the QW number at small N and 
saturates at N — > oo. The duration Tn has a minimum 
at N ~ y^ujo/To, where the transition from superradiant 
to photonic-crystalline regime occurs. The dependence 
of Vn and Tn on the indices (h,h') is determined by two 
competing effects. First, the structure factor \ fhh' \ grows 
with (h,h') [see Eq. ([%2"|)]. which enhances the structure 
reflectivity Second, the average period d increases with 
(h, h') [Eq. (|2"0|) ] leading to the increase in the flight time 
and, by that, to the suppression of the reflectivity. At 
large numbers of (h, h'), when the structure factor is sat- 
urated and close to unity, the second effect dominates. 
As a result, the energy of the reflected pulse Vn becomes 
smaller while the pulse duration Tn becomes longer. This 
explains the difference between the curves in Fig. [3] 



B. Gradient structures 

In addition to periodic or quasiperiodic structures, 
which may be characterized by the mean period, one can 
consider spatially inhomogeneous photonic systems. As 
an example of such systems, we analyze here the gradient 
structures, where the distance between the neighboring 
QWs gradually changes from the structure forefront to 
its back edge. We assume that the positions of QWs are 
given by 



Zj = {j a j)d, 



(26) 



where Q is a dimensional parameter describing the struc- 
ture inhomogencity and d is the length given by the Bragg 
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FIG. 6. (Color online) The envelopes of reflected pulses cal- 
culated for the gradient structures with various numbers of 
QWs, C = 5, and ojo/Yq = 4 x 10 4 . Figures (a) and (b) corre- 
spond to weak, T = To, and strong, T = 20ro, non-radiative 
decay, respectively. Insets show the reflection spectra of the 
same QW structures in the frequency domain. 



reflection condition (|5J). The front part of the structure is 
substantially detuned from the Bragg condition while at 
j 3> C the interwell distance tends to d and the structure 
transforms to the periodic Bragg structure. The optical 
response of the gradient structure is determined by the 
parameter £, the QW number TV, as well as the radiative 
and non-radiative decay rates Tq and T, respectively. 

The envelopes of reflected pulses pw(i) in response to 
the ultrashort excitation pulse are presented in Figs.[12a) 
andH2b) for different gradient structures. The curves are 
calculated numerically by the transfer-matrix technique. 
Similarly to the optical response of the Fibonacci struc- 
tures, the pulse envelope exhibits beatings even for struc- 
tures with rather small QW numbers N. However, the 
pattern of the beatings is different from that in the Fi- 
bonacci structures, and the pulse envelopes appear to be 
highly sensitive to the non-radiative decay rate T. The 
reflection spectra of the same gradient structures in the 
frequency domain r^v(w) are shown in insets in Figs.^a) 
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and[6[b). The spectra are much more asymmetric around 
the exciton frequency lo as compared to the spectra of 
the resonance periodic or Fibonacci structures, see insets 
in Figs. [5] and [5] The spectrum asymmetry is caused by 
the detuning of the structure front part from the Bragg 
condition. Note that the similar effect occurs in periodic 
structures detuned from the Bragg resonance ; 16 ! 38 how- 
ever, the integral reflectivity of the gradient structures 
can be larger than that of non-Bragg periodic ones be- 
cause the main part of the gradient structures is tuned 
to the Bragg condition. Another feature of the gradient 
QW structures is the narrow dip in the reflection spec- 
trum near the exciton resonance. Our calculations indi- 
cate that the dip is highly sensitive to the non-radiative 
decay and disappears at T = 0. 



V. POLARIZATION CONVERSION 

The resonant photonic crystals and quasicrystals can 
change the polarization state of an optical pulse if the ex- 
citon levels forming the photon band structure are split. 
The splitting can be induced by an external magnetic or 
electric field, strain, or intrinsic in-plane anisotropy of 
the QWs, e.g., in low-symmetry structures^ In this sec- 
tion, we consider the polarization conversion caused by 
the magnetic field applied along the structure axis (Fara- 
day geometry). The field splits the exciton level in each 
QW into two states interacting with the right- and left- 
handed circular polarized radiation. The corresponding 
resonance frequencies are given by Wq = loq ± 0c /2, 
where Ql is the Larmor frequency corresponding to the 
Zeeman splitting. The level splitting results in the dif- 
ference of the reflection spectra for the right- and left- 
handed circularly polarized light, which leads, in turn, to 
the change in the reflected pulse polarization with respect 
to the excitation pulse. For definitencss, we focus on the 
effect of the magnetic field on the resonance frequencies 
and neglect a possible change in the exciton oscillator 
strength. 

To describe the reflected pulse polarization, we intro- 
duce the Stokes parameters Si, S2, and S3, which give 
the degree of linear polarization along the axes x and 
y perpendicular to the growth direction z, the degree of 
linear polarization along the axes rotated by 45°, and the 
degree of circular polarization, respectively* 3 ^ by 

Si(t) = (\E R , x (t)\ 2 -\E R , y (t)\ 2 ), 

S 2 (t) = (E RiX (t)E* R Jt) + E R , y (t)E* RiX (t)) , 

S 3 (t) = -i(E R/x (t)E R Jt) - E R , y {t)E RtX {t)) . (27) 

Here, the angular brackets denote the averaging over the 
time domain longer than 2tt/ujo but shorter than 1/A and 
1/Oi. Note, that the sign of S3 takes into account the 
fact that the reflected wave moves in the —z direction. 
The intensity of the reflected wave is given by 



Straightforward calculation shows that the Stokes pa- 
rameters and intensity of the reflected wave have the form 

Sl(t) = RC[P + P*_](J XX - Jyy) +lm[p + p*_](J X y + Jy X ) , 



5*2(0 = Re[P+P-}(Jxy + Jyx) - lva[p + p*J\{Jxx ~ Jyy) 



I R (t) 



\P+\ 2 +\P-\\ j +J U^-IP-IV _J J 
2 \"xx~' J yyJ cy L \ u yx • J xy) j 

(29) 

where p + = p^_ + {t) and p- = /Ojv,— (t) are the response 
functions for the right- and left-handed circularly polar- 
ized radiation, respectively, and J a p {a, ft = x,y) is the 
tensor determined by the incident wave, 



c+00 



J, 



r+00 

E , a (t)e iUot dt / E*jt{t')c- iuiat ' dt' 



(30) 

In the particular case of the incident wave linearly polar- 
ized along the x axis, the Stokes parameters and intensity 
of the reflected wave assume the form 

Si(t) = Re[p_p+] J xx , S 2 {t) = Im|/»_/>* ] J xx , (31) 



The excitation of a periodic QW structure by the lin- 
early polarized pulse results in the linearly polarized re- 
sponse, whose polarization plane slowly rotates in time 
due to the Zeeman splitting of exciton states. The cir- 
cular polarization does not occur because, for periodic 
structures, |p+(t)| = |p-(t)| which follows from the rela- 
tion 



r+(w) = r*_ (2wq — u>) 



(32) 



The latter can be readily seen from the transfer matri- 
ces through the QW structures detuned from the Bragg 
condition by the magnetic field. Thus, the polarization 
state of the reflected signal can be completely described 
by the slowly varying phase tp(t), 



Si(t) + iS 2 (t) = I R (t) exp[i(p(t)] , 



(33) 



I R (t) = (\E RtX (t)f + \E Ri y(t)\ 2 } . 



(28) 



which determines the light polarization plane. Fig- 
ures [TJ^a) and[7{b) show the dependences (p(t) calculated 
numerically for structures with various QW numbers, 
weak [Fig. EJa)] and strong [Fig. [7£b)] non-radiative de- 
cays. For comparison, Fig. [TJc) presents the time depen- 
dences of intensities I R (t) of the same reflected pulses. 
One can see that, in contrast to naive expectation that 
the polarization plane rotates with the Larmor frequency, 
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FIG. 8. (Color online) (a) Time dependence of the normal- 
ized Stokes parameters determining the reflected pulse po- 
larization for the Fibonacci structure with N — 50 QWs, 
wo/To = 4 x 10 4 , T = r , and Q, L = 50r . (b) Time de- 
pendence of the reflected pulse intensity. The inset shows the 
reflection spectra for the left-handed and right-handed circu- 
larly polarized light in the frequency domain. 



FIG. 7. (Color online) (a) and (b) Time dependences of the 
phase (p, which determines the reflected pulse polarization 
plane, for the periodic structures with various numbers of 
QWs, u /T = 4 x 10 4 and U L = 50r . Solid curves are 
the result of numerical calculation after Eqs. ()31|) . dashed 
curves show the linear dependence <p{t) = Q.L.t/2. Figures (a) 
and (b) correspond to weak, Y — To, and strong, T — 20Fo, 
non-radiative decay, respectively, (c) Time dependence of the 
reflected pulse intensities. 



the dependence ip(t) is different for structures with differ- 
ent QW numbers N and can be non-linear. The simple 
law ip(t) = Q^t is valid only for structures with small QW 
numbers, TV <C yWrj/To, cjq/^l, where the reflection co- 
efficients r±(cu) are given by Lorentzian Eq. (|10j) with ujq 
being replaced by luq ± respectively. The reflec- 

tion coefficients r±(uj) of the semi-infinite structures in 
the magnetic field are described by Eq. (fT2|) . where w = 
(w-w )/A and w = (u) - uj T ^l/2 + ir)/A. At T = 0, 
they fulfill the equation r + (ui + fii/4) = r_(w — 57^/4). 
Consequently, the frequency of the plane rotation be- 
comes twice smaller and given by $1^/2. The linear func- 
tion tp(t) = shown by dashed curve describes well 
the numerically calculated dependence ip(t) for the semi- 
infinite QW structure with weak non-radiative decay, see 
Fig. 0a). The points where the phase abruptly changes 
are close to the points of minima in the pulse intensity. 
The reduction of the plane rotation frequency in the semi- 
infinite structure by the factor of two as compared to 
that in short structures can be interpreted by the for- 
mation of polariton (coupled exciton and photon) modes 
in the semi-infinite structure. The magnetic field affects 
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FIG. 9. (Color online) (a) Time dependence of the normal- 
ized Stokes parameters determining the reflected pulse polar- 
ization for the gradient structure with TV = 50 QWs, £ = 5, 
uj /r = 4 x 10 4 , r = To, and n L = 50r . (c) Time de- 
pendence of the reflected pulse intensity. The inset shows the 
reflection spectra for the left-handed and right-handed circu- 
larly polarized light in the frequency domain. 



the exciton part of the polariton only, which leads to the 
Zeeman splitting of polaritons that is twice as small as 
the Zeeman splitting of free excitons. 

Figures [5] and [5] show the time dependence of the 
normalized Stokes parameters Sj(t) — Sj(t)/In(i) (j = 



10 



1,2,3) and intensity In(t) of the optical pulse reflected 
from the Fibonacci and gradient QW structures, respec- 
tively, in the response to excitation with the short pulse 
linearly polarized along the x axis. In aperiodic systems, 
the relation (|33|) , valid for periodic structures, does not 
hold anymore. Therefore, the conversion of linear to cir- 
cular polarization becomes possible, i.e., s^it) ^ 0, see 
Figs. [HJa) and[^a). However, the efficiency of such a 
process drastically depends on the QW spatial arrange- 
ment. In particular, the conversion efficiency is rather 
low in the Fibonacci structures. This is caused by the 
fact that Eq. ([3"3"| still holds for the smooth part of the 
reflection spectra, see inset in Fig. |8jb) . Thus, the linear- 
to-circular polarization conversion occurs due to the fine 
structure of the reflection spectrum central part and can 
be theoretically described only beyond the two-wave ap- 
proximation. In the gradient structures, the reflection 
spectra are strongly asymmetric [Fig. |5Jb)], which leads 
to a more efficient polarization conversion. 



VI. SUMMARY 

We have developed the microscopic theory of the op- 
tical response of one-dimensional photonic structures 
formed by periodic, Fibonacci, or gradient sequences of 
quantum wells to ultrashort polarized optical pulses. It 
is shown that the response depends on the structure op- 
tical length as well as the QW spatial arrangement. The 
excitation of short periodic QW structures tuned to the 
Bragg condition with an ultrashort pulse leads to forma- 
tion of the reflected signal with a sharp front followed 
by a slow exponential decay. Both the signal amplitude 
and the decay rate increase with the QW number. In con- 
trast, the reflected pulse envelope in long and, particular, 
semi-infinite structures exhibits Bessel oscillations with 
the period determined by the stop band width. Moreover, 
the trailing edge of the pulse reflected from the semi- 
infinite structure has a power asymptotics, rather than 



exponential one, even in the presence of non-radiative 
exciton decay in QWs. The qualitative change in the op- 
tical response of short and long structures is caused by 
the transition from superradiant to photonic-crystalline 
regime. In the Fibonacci or gradient QW structures, the 
reflected pulses exhibit oscillations even for structures 
with rather small QW numbers operating in the super- 
radiant regime. 

The duration of reflected pulses is shown to non- 
monotonically depend on the QW number, both for 
periodic and aperiodic structures. The duration de- 
creases with the QW number in short structures, reaches 
the minimum in the area of the superradiant-photonic- 
crystallinc regime transition, and then again increases 
in long structures. This behavior suggests the optimal 
number of quantum wells for designing the fast-reflecting 
Bragg structures. 

We also have studied the conversion of light polar- 
ization in the external magnetic field applied along the 
growth direction. In this geometry, the excitation of a 
periodic QW structure by the linearly polarized pulse re- 
sults in the linearly polarized response, whose polariza- 
tion plane slowly varies in time (time-resolved magneto- 
optical Kerr effect). The plane rotation frequency in 
short structures is given by the Larmor frequency of ex- 
citons while in long structures the rotation frequency be- 
comes twice smaller, which is related to the formation of 
polariton modes. In the Fibonacci and gradient struc- 
tures, the reflection of linearly polarized pulse leads to 
both the rotation of the polarization plane and the ap- 
pearance of partial circular polarization. The linear-to- 
circular polarization conversion drastically depends on 
the QW spatial arrangement vanishing for the periodic 
Bragg structures and is more pronounced in the gradient 
QW structures. 
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